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Abstract 

We show that the parameter space of axion-hke particles can be severly constrained 
using high-precision measurements of quasar polarisations. Robust limits are derived 
from the measured bounds on optical circular polarisation and from the distribution of 
linear polarisations of quasars. As an outlook, this technique can be improved by the 
observation of objects located behind clusters of galaxies, using upcoming space-borne 
X-ray polarimeters. 

1 Introduction 

From the theoretical point of view, an observation of "axion-like particles" (ALPs) would 
provide us with valuable information about theories beyond the standard model. Indeed, 
these spinless particles — with properties similar to those of the QCD axion, but with dif- 
ferent masses and couplings — arise as testable predictions, as they are associated with the 
spontaneous symmetry breaking of new U(l) symmetries [1]. 

From the phenomenological point of view, such particles are also extremely appealing. 
First of all, depending on their masses and couplings, they can constitute a natural cold-dark- 
matter candidate; see, e.g. Ref. [2]. Moreover, as they can have a coupling with light similar 
to that of neutral pions, the existence of ALPs could help understand various astrophysical 
observations [3]. It is quite striking that these astrophysical hints seem to point towards the 
same kind of ALP, with a mass m < 10"^'' eV and a coupling to photons g ~ 10~^^ GeV~^, 
defining a whole new window of interest in the (m, g) parameter space. 

To date, the most stringent constraints on ALPs are coming from astrophysics [4], and 
are mostly related to stellar dynamics: the existence of such particles would imply the 
opening of new channels to evacuate energy, which should have an impact on stellar evolution. 
A few years ago, the limit derived from the observation of the Horizontal Branch in the 
Hertzsprung-Russell diagram for globular clusters [5] has been improved by CAST (the 
CERN Axion Solar Telescope), which tries to convert axions coming from the Sun [6] back 
into (X-Ray) photons. Discussions of the energy loss in the SN1987A supernova have also 
lead to an even more severe upper bound on their couplings to photons [7], which applies 
for nearly massless ALPs. 

The existence of such particles can also be probed using another consequence of their 
electromagnetic coupling. Until recently, it was hoped that a puzzling observation, involving 
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quasars, could also be explained by low-mass ALPs. Indeed, in some huge (~ 1 Gpc) 
regions in the sky, redshift-dependent correlations in the distribution of polarisation position 
angles for visible light from cosmologically distant quasars have been observed [8-11]. As 
one of the most interesting attributes of the mixing of axion-like particles with light is the 
induced modification of its linear and circular polarisations when it propagates in external 
magnetic fields [12,13], the presence of ALPs could have elegantly solved the problem [14-18]. 
Given the present polarisation data [19], this scenario is now strongly disfavoured, already 
at the qualitative level [20]: despite the use of a wave-packet treatment, the inclusion of 
fluctuations, and of more complex magnetic field morphologies, there is either no alignment, 
too much circular polarisation, or too much linear polarisation. Note that a recent claim of 
similar coherent orientations of polarisation in radiowaves [21] would further disfavour an 
ALP mechanism [22]. Henceforth, we shall not discuss the coherent alignment effect further. 

What we do in the present work is to use such polarisation data to constrain the properties 
of low-mass ALPs. The change of polarisation induced by the mixing is indeed a very specific 
prediction, and this is especially true when it comes to circular polarisation. Depending on 
mass and coupling of the ALPs, the mixing with light in external magnetic fields can be 
very efficient and contradict observations, as there is indeed little room for a modification of 
polarisation. Similar polarimetric constraints on these particles, based on linear polarisation, 
have already been discussed in the literature in various contexts, e.g. [14,23-25]. Being as 
conservative as possible and using present quasar polarisation data in visible light, we show 
that one can derive new constraints on ALPs using supercluster magnetic fields, and narrow 
down the parameter-space region of astrophysical interest. 



2 Polarisation of light and mixing with ALPs 

The polarisation state of any light beam Sr{z,t) = £r^{z,t)ei + £r2{,z,t)e2 — where (61,62) 
define the plane transverse to the propagation direction — is fully described by its Stokes 
parameters (which are averages over the exposure time). These can be written as 











U{z) 


= {U{z,t)) 


V{z) 


= (v(^,t)) 



(1) 



and, by combining them, one introduces two intrinsic physical quantities: the degree of linear 
polarisation and the degree of circular polarisation, respectively 



Before we use polarisation data to derive new constraints on ALPs, let us recall [20] that, 
whenever a mixing with axion-like particles takes place inside external magnetic fields, the 
polarisation of light changes. Even for initially unpolarised light, this mixing will lead to 
a modification of the linear polarisation (due to dichroism), while, for partially polarised 
light, some interplay between the circular and the linear polarisations can also occur (due 
to birefringence). These properties emerge from the fact that, given the direction of the 
transverse external magnetic field, ALPs only couple to a specific direction of polarisation 
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of light, in a slightly different^ way depending on the scalar or pseudoscalar character of the 
ALP. In this paper, we are going to particularise to the case of pseudoscalar ALPs but bear 
in mind that similar results hold for scalar ALPs. 

In regions of constant magnetic field, the evolution of the polarisation depends only on 
two dimensionless parameters [20]. One of them is the mixing angle, which unambiguously 
gives the maximum amount of polarisation that can be achieved, 

and which involves the coupling constant between photons and pseudoscalars the trans- 
verse external magnetic field strength the frequency of the light beam the ALP mass m, 
and the plasma frequency Wp. The second one determines the wavelength of the oscillations 
of the Stokes parameters, and is the product 

= , 4 

UJ UJ 

with AyU^ the difference of the squared masses of the two new eigenstates of propagation, 
and L the distance travelled. Note that this is very similar to what one obtains for neutrino 
oscillations. 



Things are more complicated in more general magnetic fields, mostly because variations 
of the electron density and of the magnetic field (both in strength and direction) imply 
variations of the above dimensionless quantities. Nonetheless, the coupling constant and the 
magnetic field strength always appear together as a product in the equations. 

In order to tackle this problem, we shall discretise the continuous magnetic field and the 
electron density in a succession of zones, and work at the amplitude level for the electromag- 
netic and axion fields, ensuring continuity at each boundary.^ The Stokes parameters are 
computed at the end of the last zone. 

Finally, note that, due to the mixing with axions, even an initially unpolarised light beam 
will be slightly linearly polarised at the end of the first zone, so that one generally expects 
some amount of circular polarisation to be created in such a field, as a result of birefringence. 



3 Polarisation data 

3.1 Properties of the full sample 

An all-sky sample of 355 good-quality measurements of the linear polarisation of visible 
light from quasars^ has been defined in a series of papers [8-10], based on data obtained 
in Refs. [26-30]. In that sample, only objects with a polarisation degree > 0.6% and 
Piin/o"piin > 2 were considered. In the present work, we extend the sample to include low- 
polarisation objects with < 0.6% measured with uncertainties ap^.^_^ < 0.3%. Since pun is 



^Scalars only couple to photons with a polarisation perpendicular to the magnetic field, while pseu- 
doscalars only couple to parallel ones. 

^Note that we neglect reflected waves, which have an amplitude of order (y^^^- 

^Here, "quasar" refers to "high-luminosity AGN." 
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Figure 1: Comparison of the distributions of debiased linear polarisation between non-BAL 
quasars {left) and BAL quasars [right), for objects taken in the direction of region Al. 
Clearly, BAL quasars are often more linearly polarised than quasars without broad absorption 
lines. 



positive definite, it is biased at low signal-to-noise ratio. A reasonably good estimator of the 
true (debiased) polarisation degree is computed using = {p\^ — when > cjpjj^, 

and piin = when < cr^jj^ [31]. In the rest of this paper, we only consider when we 
discuss the linear polarisation. 

Circular polarisation measurements using a Bessell V-filter^ [32] for some objects in the 
sample are available in [19], together with a compilation of previous measurements [26,33-35], 
most of them in white light (with no filter^). Out of 21 V-filter measurements reported 
in [19], all but two are compatible with no circular polarisation at 3(T. The two positive 
detections concern highly linearly polarised blazars (both with pw^ > 22%), which could be 
intrinsically circularly polarised [19], as in the case of a highly linearly polarised BL Lac 
object (piin > 26%), for which a non-zero circular polarisation in white light has also been 
previously observed [19,35]. These special cases set aside, light is essentially not circularly 
polarised, be it white or V-filtered. 

3.2 Subsample and criteria used to obtain constraints 

Since our goal is to derive constraints from the polarisation induced by axion-like particles, 
we restrict our sample to spectroscopically defined classes of quasars known to have the 
smallest intrinsic polarisations. 

We discard the following objects from the sample of quasars with linear polarisation 
measurements: known radio-loud quasars [i.e., classified as such in the NASA/IPAC Extra- 
galactic Database (NED) or objects from references [30]), near-infrared (2MASS) selected 
objects, in particular those from [29], and Broad- Absorption-Line (BAL) quasars (classified 
as such in NED). A comparison of non-BAL measurements to the polarisation distribution 

■^This filter is centered around A ~ 547.6 nm, with a full width at half maximum of 113.2 nm. 
^For instance, the photomultipliers used in [33] have a broad spectral response, from 185 to 930 nm. Note 
however that there is an atmospheric cutoff for wavelengths below 330 nm [36]. 
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Figure 2: Distribution of the degree of circular polarisation for the central values of the 16 
quasars in the direction of region Al, measured in the V-filter. 



of BAL quasars is shown in Fig. 1, demonstrating that differences are significant; they are 
also preserved by any mechanism possibly affecting light on the line-of-sight, see Ref. [9]. 

In the sample of circularly polarised quasars, we discard the two spectroscopically iden- 
tified BL Lac objects, more prone to be intrinsically polarised. 

The reduced sample then mostly consists of radio-quiet / optically selected non-BAL 
quasars. Such objects are known to be linearly polarised at most at the 1% level [26,27]. 

As the observation of null circular polarisation is a much more stringent constraint for 
axion-like particles in V-filter than in white light [20], we focus on quasars located in the 
direction in which most of these V-filter data were taken, i. e. in the direction of the so-called 
"region Al" [8-10], i.e. towards the center of the (local) Virgo supercluster. Keeping only 
objects with right ascensions between 168° and 218°, we are left with a final sample of 55 
quasars with measured linear polarisation in white light or in V-filter and of 16 quasars with 
measured circular polarisation in V-filter (all compatible with zero at 3a). The distribu- 
tion of the circular polarisation degree associated with the central values of the 16 circular 
polarisation measurements is shown in Fig. 2. 

We know that these quasars intrinsically emit polarised light: otherwise we would simply 
not see a difference in the distribution of degrees of linear polarisation depending on the 
spectroscopic type. However, we cannot access the initial distribution of polarisation. For 
this reason, in order to avoid any overestimation of the final polarisation generated by the 
mixing, we work with initially unpolarised light. We conservatively allow the observed 
linear polarisation of the quasars to be only due to the interaction with ALPs. Under this 
hypothesis, one can impose that the linear polarisation generated by the ALP-photon mixing 
model does not exceed the observed one. 

The idea is therefore to calculate the polarisation predicted by the mixing with ALPs 
inside external magnetic fields encountered on the way to Earth and to compare it with data. 
For both the linear and the circular polarisations, we thus compute the probability that the 
polarisation (p**^) due to ALPs is smaller than the observed one {p°^^) as: 

P = Ar(p°bs > pti^)/N,,,,^^ (5) 
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by taking the ratio of the number of measurements in bins with a polarisation higher than 
p**^, to the total number of measurements A^totai- 

4 Model for supercluster magnetic fields 

4.1 Structure 

On the observational side, what is typically assumed in the literature for magnetic fields inside 
clusters or superclusters is a cell-like structure, see e.g. [37-39], possibly with an additional 
weaker background field coherent over several cells. One usually considers a collection of cells 
of a given size, and lets the magnetic field direction change from cell to cell while keeping 
the same field strength |i?ceii|,^ which averages to a smaller value at the supercluster scale. 

In such a configuration, one sees different fields when looking in different directions. 
Because of this, and because the exact magnetic structure is not known inside the Virgo 
supercluster, we assume that the effect of these irregularities is equivalent to considering 
that light coming from different objects essentially passes through random fields. 

Additionally, we perform general tridimensional rotations of the magnetic field direction 
from cell to cell; it will thus pick up a longitudinal component most of the time. As mentioned 
in Sec. 2, only the projection of the magnetic field onto the tranverse plane is relevant for 
ALP-photon mixing. Hence, we allow values of the transverse magnetic field strength B 
much smaller than |-Bceii|- 

We also consider fluctuations of the electron density from cell to cell and of the size of 
each cell. More precisely, we let these quantities fluctuate by 50% up or down. Note that, 
as we allow for fluctuations of cell sizes, only the total distance Ztot is well-defined from one 
magnetic field realisation to another, not the number of cells, as we stop adding cells as soon 
as we reach Ztot- 

Finally, note that the size of the supercluster (~ 20 Mpc) is small enough so that cosmo- 
logical redshift can be neglected. As the mixing is more efficient at higher frequencies (see 
Eq. (3)), redshift effects could only increase the polarisation produced in any case relevant 
to this study and make our bounds stronger. Moreover, as we only take into account the 
influence of the last magnetic field region, we again underestimate the amount of polarisation 
due to ALP-photon mixing. We thus stress that our constraints are conservative. 

4.2 Field strength 

Magnetic fields in galaxy clusters are well established [38,41,42] and have typical strengths 
of several fiG, although they can be as large as 40 fiG. Our knowledge of these fields is 
sometimes refined enough that not only the field strength and the coherence lengths are 
known but also hints about the orientation are available [43]. In our case, even though 
region Al is roughly centered on the Virgo cluster, most of the quasars are seen through a 
larger domain, and we need to model the magnetic field at the scale of the supercluster. 

There is evidence for a sizeable magnetic field in the surroundings of galaxy clusters [44], 
and some results are also available for superclusters. For instance, it has been inferred from 

^Note that this kind of model is also what is usually used when one discusses magnetic fields beyond 
supercluster scales; see for instance [40]. Such faint cosmological magnetic fields have also been used to 
obtain limits on ALPs [25]. 



6 



rotation measures in radio wavelengths [37,38] that the structure of the magnetic field within 
the local supercluster plane (centered on the Virgo cluster, the direction we focus on) can 
be described as a collection of ~2 nG magnetic field zones coherent over ~100 kpc. If one 
interprets these data with a magnetic field coherent over the supercluster scale, or assumes 
much larger coherence lengths, one then obtains field strengths about 5 times smaller. A 
similar larger-scale magnetic field is also considered in Ref. [39], which reports upper values 
for typical magnetic field strengths in the Hercules and in the Perseus-Pisces superclusters 
of the order of 0.3 ± 0.1 /iG (over 800 kpc) or 0.4 ± 0.2 /iG (over 400 kpc). 

5 Minimal constraints 
5.1 Method 

Our constraints on ALPs are derived iteratively for different points in the (m, g) parameter 
space: a given couple of parameters corresponds to a given hypothetical axion-like particle. 
For fixed ALP parameters, we generate random configurations, defined as the set of trans- 
verse field strengths and directions, of cell sizes, and of electron densities. For a given value 
of the frequency u, we have the solutions for the two electric fields of the radiation after a 
propagation in each "patch" of such a configuration [20]. We use initially unpolarised light 
beams (as discussed in Sec. 3.2), calculate the four Stokes parameters at the end of the last 
magnetic field cell, and compare the final polarisation with data. 

Now, as explained in [20], the expected circular polarisation due to ALP-photon mixing 
can be reduced when one considers wave packets, or averages over frequencies. This effect, 
though typically small in V-filter, is taken into account when deriving our limits. 

To implement this, for a single light beam, we completely generate a random configura- 
tion. We then integrate over u the monochromatic results for the Stokes parameters at the 
end of the last zone of magnetic field, with the spectral response of the V-filter as a weight- 
ing function. This corresponds to one iteration: it gives the values of the Stokes parameters 
through the V-filter for a given source at the end of one configuration. 

After that, we compare the generated linear and circular polarisations with the data 
histograms^, associate to each of them a probability given by Eq. (5), respectively Pun and 
Pcirc, and obtain the final probability V to be compatible with data for this particular con- 
figuration by multiplying these two individual probabilities. We then repeat the procedure 
with a new random configuration. 

We do this many times (5000 times) to minimise the infiuence of statistical fiuctuations, 
and average the value of V over the configurations as we do not have access to the actual 
configuration for each measurement. We thus integrate it out and give the average probability 
for a given ALP (m, g) not to exceed the observed polarisation in a random configuration. 
Then we start over for a new couple (m,g). 

Finally, we summarise this information by saying that parameters are "excluded at la" , 
when the average probability for this particle to produce too much polarisation compared 
with data is 68.3%; and similarly for 2a (95.5%) and 3a (99.7%). In Appendix A, we show 
two theoretical distributions of polarisation, associated with points excluded at 2a and 3cr; 
let us stress that this is only for illustration, as the technique described above to obtain 
constraints does not directly involve such distributions. 

''The results are not affected by the presence of white hght data for hnear polarisation. 
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5.2 Results 
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Figure 3: Exclusion plot. The total distance considered here is Ztot = 10 Mpc; the average cell 
size is 100 kpc; the norm of the total magnetic field in each cell is |-Bceii,o| = 2 fiG; the average 
electron density is = 10~^ cm~^ (corresponding to a plasma frequency Up = 3.7x IQ-^^ eV). 



The results we obtain are shown in Fig. 3. The upper limits obtained by CAST and from 
the energy-loss considerations associated with SN1987A are also shown, and confirmed for 
nearly massless ALPs. 

Let us first stress a very important point: as there is only one physical magnetic field 
scale in the problem, namely the total field strength inside a cell |-Bceii|; this exclusion plot 
remains the same for any other value of the magnetic field inside the cells. Indeed, as 
stressed in Sec. 2, only the product of the magnetic field with the coupling constant g enters 
the equations. Therefore, while there are uncertainties about the magnetic field strength, 
as discussed in Sec. 4.2, the exclusion limit on g can be rescaled once our knowledge of the 
magnetic field inside the supercluster improves. 

In Fig. 3, we observe two features: a dip and a plateau for small ALP masses. Both 
can be understood from the expression of the mixing angle, Eq. (3). The position of the 
dip is determined by Up] for ALP masses very close to the plasma frequency, the mixing 
is maximal, and so is its effect on polarisation. As for the plateau, for ALP masses much 
smaller than the plasma frequency, the mixing angle is essentially independent of m: 

(m«<a;p) 1 / 2gBuj\ 

fmix = -atan ; (6) 



2 V -00 



8 



this exclusion thus holds even for massless ALPs. Now, the fact that higher values of the 
coupling g are more constrained is natural, as the mixing is more efficient and, thus, more 
polarisation is produced. 

Here, for the average electron density inside the supercluster, we have used Ue = 10~® cm~'^ 
(which gives a plasma frequency^ Up = 3.7 x 10~^^ eV), as in previous works related to 
ALP-photon mixing in similar conditions, e.g. [15,45]. Nonetheless, the properties of the 
intrasupercluster medium are not well-known, including the value of the average electron 
density. Searches for gas in the Shapley supercluster [46] have bounded rie to be less than 

5 X 10~^ cm~^, and a subsequent tentative detection, for the local supercluster [47], gave 
TT-e ~ 2.5 X 10~^ cm~'^. More recent constrained N-body simulations of the local supercluster 
in the direction we are looking at (at galactic latitudes b > 30° and distances of about 10 Mpc) 
lead to values for the plasma frequency between Up = 1.7 x 10"^"^ eV and Up = 5.5 x lO"^'' eV 
for matter overdensities between 5g = and 6g = 10 [48]. 

Now, if we take an extreme case and allow the intrasupercluster electron density to be 
one order of magnitude larger on average than what we have considered in Fig. 3 {i.e. allow 
fluctuations, see Sec. 4.1, up to values as large as rie = 15 x 10"^ cm~'^), we have checked 
that, while our limits would then change, they would still improve the current bounds on 
ALPs. For smaller values of the electron density Uc < 10~^ cm~^, our limits are stable, and 
would in fact be slightly more stringent. 

We have checked the stability of our constraints under the change of the total magnetic 
field size Ztot- For total magnetic field sizes of 5 Mpc and 20 Mpc, the exclusion plot obtained 
essentially does not change its shape compared to Fig. 3, but is shifted along the y-axis. 
For instance, the 2(j-limit we obtain for nearly massless ALPs is < 2.5 x 10~^^ GeV~^ for 
ztot = 10 Mpc, and becomes ^ < 3.2 x 10"^^ GeV"^ and ^ < 2 x lO'^^ GeV-^ for ^tot = 5 Mpc 
and Ztot = 20 Mpc respectively. 

We have also checked that our constraints are stable: a) under the addition of a uniform 
background field (typically weaker: |-Buniform| ~ 0.4 /iG), which would allow some correlation 
between cells over the supercluster scale;^ b) if we use the histograms of linear polarisation 
for non-BAL quasars measured in white light and V-filter, or exclusively in V-filter to define 
our probability; c) if we include the measurements of circular polarisation for the two BL 
Lac objects; d) if we do not allow for fiuctuations of the size of the cells and of the electron 
density from cell to cell. 

6 Conclusions and outlooks 

This paper extends previous proposals to use polarisation data to constrain ALP parame- 
ters [14, 23] by using a realistic magnetic field and an average over frequencies. We have 
presented constraints derived simultaneously from linear and circular polarisation measure- 
ments, sticking to a conservative approach where we allowed the observed polarisation of 
non-BAL quasars to be only due to ALPs. Using reported properties of the intrasuperclus- 
ter medium, we have showed that current bounds on nearly massless ALPs are improved. 
Indeed, for a large portion of the parameter space region of astrophysical interest, the po- 
larisation produced by the mixing with such particles in this medium would be too large. 

^The plasma frequency is given by uip ~ 3.7 x 10^^'*-\/?ic/10~^ cm"^ eV. 

^Then, of course, one cannot rescale the limits derived for |-Bccii| and |i?uniform| changed independently. 
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While our constraints can be straightforwardly rescaled when our knowledge of the mag- 
netic field strength improves, the weak point of the method is of course the uncertainties on 
the electron density and on the field strength, and more reliable constraints could be obtained 
using galaxy clusters where these properties are better determined. As the plasma frequency 
is much higher at the cluster level, an effect in visible light would however be suppressed, 
except at the resonance. Using polarisation data in X-rays instead, one could take advantage 
of galaxy clusters to constrain ALP properties; this could be achieved with instruments such 
as the polarimeter onboard of the GEMS satellite [49], which will operate between 2 and 
10 keV. Along with other authors [50], we therefore stress that the polarisation properties 
of very-high-energy photons is a promising tool to search for ALPs. 
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A Illustration of theoretical distributions 

We show distributions for the linear and the circular polarisations obtained after 5000 dif- 
ferent magnetic field realisations. We reproduce what is obtained for two points taken from 
Fig. 3, both for an ALP mass m = 10~^° eV, but for different values of the coupling. 

In Fig. 4, we show the theoretical distribution of polarisations obtained for initially un- 
polarised light for a coupling g = 2.5 x 10"^'^ CeV~^, which is the lower bound of parameters 
excluded at 2a for nearly massless ALPs. In Fig. 5, we show the same for a coupling 
g = 5 X 10~^^ CeV~^, i.e. the 3a lower bound for nearly massless ALPs. 

In both case, we see that ALP-photon mixing produces much more polarisation than 
observed (Figs. 1 and 2); in particular, the circular polarisation is very large. For values of 
the coupling close to the current bound from SN1987A, the linear and circular polarisations 
go above 50% of polarisation in some configurations while we would expect at most around 
1% for the linear polarisation and a very tiny circular polarisation. 
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Figure 4: Theoretical distributions of polarisation for 5000 randomly generated configu- 
rations for initially unpolarised light. Here, m = 10"^'' eV and g = 2.5 x 10~^^ GeV~^; 
parameters and fluctuations are the same as in Fig. 3. 
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Figure 5: Theoretical distributions of polarisation for 5000 randomly generated configura- 
tions for initially unpolarised light. Here, m = 10~^° eV and g = 5x 10"^'^ GeV~^; parameters 
and fluctuations are the same as in Fig. 3. 
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